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Abstract 

We present the full set of evolution equations for the spatially homogeneous cosmologies 
of type Vlh filled with a tilted perfect fluid and we provide the corresponding equilibrium 
points of the resulting dynamical state space. It is found that only when the group parameter 
satisfies /i>-la self-similar solution exists. In particular we show that for h > — | there 

exists a self-similar equilibrium point provided that 7 € (^5+3^^^' l) whereas for h < — ^ 

the state parameter belongs to the interval 7 G ^1, ^^^ZS -^ . This family of new exact 
self-similar solutions belongs to the subclass = having non-zero vorticity. In both cases 
the equilibrium points have a five dimensional stable manifold and may act as future attrac- 
tors at least for the models satisfying — 0. Also we give the exact form of the self-similar 
metrics in terms of the state and group parameters. As an illustrative example we provide 
the explicit form of the corresponding self-similar radiation model (7 = |), parametrised by 
the group parameter h. Finally we show that there are no tilted self-similar models of type 
III and irrotational models of type VI^. 

KEY WORDS: Exact Solutions; Perfect Fluid Models; Self- Similarity. 



1 Introduction 

Although on a sufficiently large observational scale, the present state of the Universe is described 
by the Friedmann-Lemaitre (FL) model which is isotropic and spatially homogeneous, there are 
potential problems mainly regarding the observed local structures of our "lumpy" Universe which 
cannot be explained within the class of FL models. Therefore more general cosmological models, 
which in some dynamical sense, are "close" to FL but not isotropic in local scale, can be used in 
order to answer many open and important questions. For example it is of interest to understand 
the presence, the form and the evolution of small (local) density and expansion anisotropies in 
the Universe or to investigate the constraints that measurements of temperature anisotropies 
are able to impose on the curvature of space-time. In addition it is important to classify all 
possible asymptotic states near the cosmological initial singularity (i.e. near the Planck time) 
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and into the future that are permitted by the Einstein's Field Equations (EFE) with a view to 
explaining how the real Universe may have evolved. 

In this point of view the simplest (anisotropic) generalisation of the FL universes are the Spa- 
tially Homogeneous (SH) cosmologies admitting a G3 group of isometries acting on 3-dimensional 
spacelike hypersurfaces C. From cosmological point of view, it is of importance to study the 
evolution of vacuum or models filled with a gamma-law perfect fluid matter source having an 
energy-momentum tensor of the form: 



where fl,p = (7 — 1)/! are the energy density and the pressure measured by the observers 
comoving with fluid velocity u a {u a u a = —1). Since in SH models there is a preferred unit 
timelike congruence n a (n a n a = — 1) normal to the spatial foliations C we can divide them into 
non-tilted p0| and tilted [2] models according to whether the fluid velocity u a is parallel or not 
to the timelike direction n a . 

In the last two decades, the study of SH models is heavily based on the qualitative analysis 
of the resulting system of the (induced) first order ordinary differential equations. Using the 
so called orthonormal frame formalism (pioneered by Ellis pQ) which is based on choosing a 
frame tetrad invariant under the group of isometries (thus ensuring the spatial independency 
of the kinematical and dynamical quantities of the models) and a set of expansion-normalized 
variables, the evolution and constraint equations, followed from the EFE, become an autonomous 
system of decoupled first order differential equations which can be studied with the aid of the 
well-established theory of dynamical systems j3]. The evolution of a specific model is studied in 
the so called dynamical state space which represents the set of all the physical states (at some 
instant of time) of the corresponding model 0. Under this perspective of studying SH models, 
equilibrium points (i.e. fixed points) and their stability, of the resulting dynamical system, play 
an important role in the description of the asymptotic behavior (into the past/future as well as 
in the intermediate times of their evolution) since they may represent past or future attractors 
for more general models. 

Vacuum and non-tilted perfect fluid SH models have been extensively studied in the liter- 
ature |S1 EJ 03 E3 13 |H1 El revealing new and important features of the SH models like e.g. the 
asymptotically self- similarity breaking and the divergence of the Weyl curvature at late times 
for type VIIo models, providing a solid counterexample to the isotropisation conjecture (i.e. 
shear isotropisation implies a corresponding result for the Weyl curvature scalar). On the other 
hand it is natural to expect that the behavior of SH models will be modified accordingly by the 
presence of a tilted fluid velocity leading also to new interesting phenomena. 

The first step of qualitative analyzing tilted models has been done for Bianchi type II models 
|10| . In particular it has been shown that 7— law tilted perfect fluid models, are future asymptotic 
to the Collins-Stewart non-tilted model ^Tj when | < 7 < ^, consequently these models do 
not isotropise and the angle of tilt becomes negligible at late times. At the value 7 = -S the 
tilt destabilise the Collins-Stewart model and there is an exchange of stability with the self- 



the asymptotically extreme tilt for models where the state parameter 7 belongs to the interval 



Recently it was shown that the self-similar equilibrium points of Bianchi type VIq models 



Tab = (A + P)u a u h + pg a b 
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play a similar role in the asymptotic behaviour of generic models. For example it was found 
|12| that at the value 7 = | the tilt destabilise the Collins solution |13| and a family of models 
satisfying = are future asymptotic to the Rosquist and Jantzen self-similar model |121 114j 
for 7 £ f|, However it has been shown that generic models (i.e. those satisfying ^ 0) 
are not asymptotically self-similar and may be extreme tilted at late times for | < 7 < 2 

m- 

In the case of class B tilted models less information is available due to the increased com- 
plexity of the evolution equations. Recently the whole family of Bianchi class B models has been 
studied and some results concerning the stability of the non-tilted equilibrium points have been 
given [Hil (see also [d EE! ) • 

Motivated from the above facts, the goal of this work is to present the full set of evolution 
equations, the vacuum, non-tilted and tilted equilibrium points for the important class of Bianchi 
type VI/,, models and interpret, for some of them, their geometric and dynamical properties. 

An outline of the paper is as follows: section 2 reviews and presents the basic results concern- 
ing the set of equations which describes the dynamics of the tilted perfect fluid SH models. By 
specialising to the case of Bianchi type VI/,, models, we provide the complete set of the evolution 
equations and we identify the resulting dynamical state space. In section 3 we find the complete 
set of equilibrium points and for the case which are represented by self-similar models, we give 
all tilted perfect fluid models admitting a proper Homothetic Vector Field (HVF). Finally in 
section 4 we summarise and discuss the implications some of the obtained results. 

Throughout the following conventions have been used: spatial frame indices are denoted 
by lower Greek letters a, /?,... = 1,2,3, lower Latin letters denote space-time indices a,b,... = 
0, 1, 2, 3 and we use geometrised units such that 8ttG = c = 1. 



2 Dynamical state space of tilted perfect fluid Bianchi type VI/j 
models 

In SH tilted perfect fluid models, the autonomous differential equation governing their evolution 
can be written in the form: 

£ = /W (2.D 

where x is the state vector representing the set of all the physical variables that describe the 
dynamics of the corresponding model, /(x) is a polynomial function of the state vector and r 
is the dimensionless time variable defined by: 

dt 1 dH , , 

where q, H are the deceleration and Hubble parameter respectively. 

In ^0] and using the orthonormal frame approach, the EFE are reformulated in terms of 
the components of the shear tensor of the normal timelike congruence n a , the spatial curvature 
of the orbits of the G3 isometry group and the spatial part of the tilted fluid velocity u a . The 
evolution equations for the type VI/i models can be found by specialising the set of equations 
given in |10| : 

E' a/3 = -{2-q) £ Q/3 + 2e»\j:^R v - S a p + U a/3 (2.3) 
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n' = nG 



-i 



A a = qA a — Ti^A^ + e^ v a A^R v 
2Gq - (3 7 - 2) - (2 - 7) - 7 £ |H ,t;'V + 2 7 ^^ 



(2.4) 
(2.5) 
(2.6) 



[l-(7-l)v 2 ] 



{(3 7 - 4) ( 1 - v 2 ) + (2 - 7) E 7 i«V + 



+ 



(2 - 7) - (7 - 1) (l - v 2 )] Apv 13 } - S Q % + 



where a prime denotes derivative w.r.t. r. 

The above system is subjected to the algebraic constraints: 

n = l-T, 2 - K 

3 7 G- 1 0^ Q = 3£ a % - eJ^E/Nfr 

where we have set 

G = 1 + (7 - 1) v 2 
and the deceleration parameter is given by the relation: 



q = 2£ 2 + 2 G ~ lf] [( 3 7 ~ 2 ) (! - ^) + 2 V 

= 2 (1 - K) - \g~ 1 Vl [3(2-7) (l -v 2 ) + 2 7 v 2 



(2.7) 



(2.8) 
(2.9) 

(2.10) 



(2.11) 



Using the freedom of a time-dependent spatial rotation, we may choose the orthonormal tetrad 
to be the eigenframe of N a p therefore the contracted form of Jacobi identities N a pA@ = 
implies: 



N, 





N 2 
N 3 



A a = Arfl. 



(2.12) 



The evolution equation of N a n can be used to express the angular velocity R a of the spatial 
frame in terms of the shear variables: 



til — TT^23) -n-2 — ^13) -K3 



->12- 



(2.13) 



iV 2 - jV 3 

In addition equations (|2.4|) and (|2.5j) have a first integral which is used to express the component 
A\ in the well known form: 

A\ = hN 2 N 3 . (2.14) 
Following XOj we introduce the shear variables: 



= - (£22 + ^33) , S- — -^j^ (£22 - £33) 



(2.15) 
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y/3' 



J 23> 



— =Eio, Ei 3 — > — = Xi 3 . 
^3 ^3 



(2.16) 



In the case of type VI^ models we have h < and N 2 N 3 < 0. With these identifications 
we obtain the following set of evolution equations for the basic expansion-normalised variables 
x = (E + ,E_,Si,S 3 ,Si3,iV2,iV3,7; a ): 



Si 



(iV 2 - iY 3 ) 2 - 18 (E? 3 + E|) 7 (2«? - t; 2 2 - ^|) 



2G 



E' 



(2 - 9 ) E_ 



+\/3 (2E 2 - E 2 3 + El) + 



iV 3 2 ) 4 v / 3A r 3 Ef 



iV 2 -iV 3 



+ 



v / 3^7( 



U 5 — US 



2G 



+ 2 V3E_ + 1 



Ei + 



N 2 -N 3 

| ^ VhN 2 N 3 (N 3 - N 2 ) + 6E 13 E 3 | ygWjr 



G 



S'u = (q + V3S_ - 3S+ - 2) E 13 



N 2 - N 3 ' G 

2v / 3-/V 2 E 3 £i y/3viv 3 flj 
N 2 -N 3 + G 

jV 2 ' = (q + 2V3E- + 2E+) iV 2 

iV 3 ' = (g - 2^2- + 2E+) iV 3 

and the evolution equation (|2.7[) for the frame components of the tilted fluid velocity. 
The algebraic constraint (|2.9j) reads: 



\/3Ei (iV 3 - iV 2 ) + e^A^s^ + = 



G 



V3E 13 A^ 3 + 3V3VhN 2 N 3 Z 3 - = 



G 



V^Es^ - 3\/3V/iiV2JV 3 Ei3 + = 0. 



We note that the shear scalar E = — s- 2 



E 2 = Ei 



- and the spatial curvature K are: 
Y?_ + S 2 + E 2 + E 2 3 



if = /iiV 2 iV 3 + 



(iv 2 - jv 3 ) 
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(2.17) 



(2.18) 



(2.19) 

(2.20) 

(2.21) 
(2.22) 
(2.23) 



(2.24) 

(2.25) 
(2.26) 

(2.27) 
(2.28) 



therefore the inequality Q > and the constraint ()2.8|) imply that the state space T> C R' is 
bounded (we recall that N 2 N 3 < 0). 
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3 Determination of the equilibrium points 



Equilibrium points of the autonomous differential equation (\2.1\i play an important role in the 
evolution of the SH models since they determine various stable and unstable invariant submani- 
folds of the state space T>. These points can be found from the solution of the algebraic equations 
/(x) = and (|2.24|) - l|2.26|l which we now list in the following subsections. We note that the 
case of type III models is included by setting (whenever is appropriate) h = — 1. 



3.1 Vacuum Equilibrium Points 

1. Kasner Circle K, HI 



N 2 = N 3 = 0, v a v a = 0, £ 2 = 1 

Si = S13 = S 3 = 0, q = 2. 



£ 2 + £ 2 = 1 



2. Kasner Line with tilt tC^ lt E 



N 2 = N 3 = 0, v a v a = vl < 1, v x = v 2 = 0, X 2 = 1 



£+ = v^S- + 37-4, £_ 



V3(3 7 -4)± ^3(2-7) (37-2) 



Si = S 3 = Si 3 = 0, q = 2, 



< 7 < 2. 



3. Kasner Circle with extreme tilt tC ex treme |10| 

AT 2 = at 3 = o, = vl = 1, «! = v 2 = 0, S 2 = 1 

S 2 +S 2 =l, Si = S 3 = Si 3 =0, q = 2, 0< 7 <2. 



^. Collins Vacuum Plane Wave Arc C{VIh) (h ^ — 5) 

-12/tiV 2 iV3 = 2^7" (^2 - N3) 2 + 3 + (iV 2 - N 3 ) 2 - 6, u a t; a = 0, 



V3 + J- {N 2 -N 3 ) 2 + 3 



■1 2 



12/iiV 2 iV3 + (N 2 - N 3 y 



12 (N 2 -N 3 y 



V3 + J- {N 2 -N 3 ) 2 + 3 



2^3hN 2 N 3 Z + 
(N 2 -N 3 ) ' 



£_ = E13 = S 3 = 0, 



VS V - (N 2 - N 3 ) 2 + 3 + ^3 
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< 7 < 2. 



5. Vacuum plane wave with tilt 1 Mt u {VI h ) (h ^ 



N 2 = N 2 , N 3 = N 3 , v a v a 



(N 2 -N 3 ) 2 + 3t V3(3 7 -5) 



12hN 2 N 3 ( 7 - 1)' 



{N 2 - N 3 ) 2 + 3 ± V3 



2 r 



12/iiY 2 ^3 + (iv 2 - N 3 y 



12 (N 2 -N 3 ) 2 



V3 



(N 2 - N 3 ) 2 + 3 - V3 



J 13 



"U2 = ^3 = 



(N 2 - N 3 ) V3 



V3 



±J- (N 2 -N 3 ) 2 + 3-V3 (37 - 5) 



6y/hN 2 N 3 (1 - 7) 



V3 



x/3± J-(N 2 -N 3 ) 2 + 3 



(N 2 - N 3 ) 4 



12N 2 N 3 



±2^/3 J- (N 2 - N 3 ) 2 + 3-(N 2 - N 3 ) 2 + 6 



We remark that the state parameter 7 is constrained via the inequality 1 — v 2 > 0. 



6. Vacuum plane wave with extreme tilt M. e xtreme{VIh) (h / — g). 

Same as the case M.f ilt (VIh)- However the state parameter can take any value in the interval 
(0,2). 

1 It appears that this form of the Collins type VI^ plane wave solution has been also given in .18 using, however, 
a different notation. 
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3.2 Non Vacuum Equilibrium Points 

1. Flat Friedmann-Lemaitre Equilibrium Point T jTUj 

N 2 = N 3 = 0, v a v a = 0, S 2 = 

S + = £_ = Si = Si 3 = S 3 = 0, q = 2 , = 1 

< 7 < 2. 

2. Collins-Stewart type II non-tilted Equilibrium Point CS(II) [2] 



2 



4 ' " ' 64 

S + = lg 37 ' S - = — -> Si = Si 3 = S 3 = 

37-2 2 « 3(6-7) 

g = — — , -<7<2, n = -i — -. 

* 2 ' 3 ' 16 
5. Hewitt type II tilted Equilibrium Point Vtut(II) jHH 



_ 0, iV 3 = 3, (7-2)(37 - 4)(5 7 - 4) ; _ (37 - 4) (7 7 - 10) 



18 - 177 ' (117 - 10) (57 - 4) ' 

V 

177 - 18 " 



E 2 _ (37 ~ 4) (9 7 2 - 20 7 + 12) s _ 9 7 - 14 ^ _ ^3(57-6) 



3(7-2) (7 7 - 10) (117 ~ 10) v v n 
Sl3 = V 16(18-177) ' Sl = S3 = ' 

_ 3 (2 - 7) (21 7 2 - 24 7 + 4) 
4(177 - 18) 



(3 7 - 4) (7 7 - 10) 3 7 - 2 10 

^ = V (ii7-io)(5 7 -4r q= — y <7<2 - 

4- Type II Line of tilted Equilibrium Points Ctut(H) JO] 



2 (27^ + 2) (17 -54b 2 ) 6 (27b 2 + 1) (27b 2 + 2) 

2 "V 57 ' 3 " ' Q " (54b2 - 17) (81b 2 - 32) ' 

2 (2 -3b 2 ) (27b 2 + 2) V(27b^ + 1) (32 - 816*) 

19 ' 3^57 



2^3 6 (27b 2 + 1) (27b 2 + 2) 

£_ = — , S + = Si = 0, w 3 ^ 



(17 -54b 2 ) (32 -81b 2 ) 



4 n 2916b 4 - 1215b 2 + 236 ,,, 2 14 

-, 0.= , b < — — , 7= — . 

3 342 'M 3^/3' ' g 



8 



5. Type II Extreme tilted Equilibrium Point Vextreme(H) jlOj 



6\/19 

N 2 = ^—,N 3 = 0, v a v a = l, 
iy 

2 _28 10^57 2^3 

S "57' S3 -"^7T' Sl3 "^"' 

S- = --g-> S + = Si = 0, U 3 = 1 

4 ^20 
?=g, n = l^> 0< 7 <2. 

Collins Type VI \ non-tilted Equilibrium Point C^{VIh) [S] 



N^ ^—'l^- 2 \ N i = -N 2 , ^„ = 0, 
tf- Or -»)'('-») , E 3 = E 13 = E_=0, 

lb 

2-3 7 ±V3(2-3 7 )^ 
S + = ^' Sl " 4 > 

37-2 3[/t(3 7 -2)- 7 + 2] 2 2(l-/i) 
9= ^-' ° = 4 ' 3- 7 -TT^W 

7. Type VI h tilted Equilibrium Point CtUt(VIh) (h = —hf and /ii ^ —1/3): 

-2/i l7 [- (3/n + 1)] 1/2 {3/i 2 (4g 7 - /?) - 4/n [g (5 7 - 4) + 2 (7 - 2)] + 3/?} 1/2 iV 3 

= \/3{6/if 7 [{q + 1) (2g 7 - 0) (6q (7 - 1) + 5 7 - 6)] - ft? [ 7 3 (l20g 3 + 16g 2 - 28q + 49) + 

+ (q + 1) (-2 7 2 (144g 2 + 32g + 23) + 60 7 (g + l)(4g + 1) - 72(g + l) 2 )] + 

+h 2 1 (3[j 2 (8q 2 - 66g - 41) + (g + 1) (4 7 (g + 18) - 12(g + 1))] + 

/ii/? 2 [2g (4 7 - 5) + 13 7 - 10] - /3 3 } 1/2 

a = 3/? 2 {/if [6g (7 - 1) + 5 7 - 6] - 2hx [q (3 7 - 2) + 7 - 2] + g} 
U Va 4/i 2 iVf 7 2 {3hj (4 97 - /?) - 4/ii [g (5 7 - 4) + 2 (7 - 2)] + 3/3} 

x/3{2/iiiV3 [C + ^ 2 7 2 (7 - 1)] + 2iV 3 C + (/? + 2 - g)} 



S 



3 



12C^ 2 7 



y y y V* frl (g + 1) (5 7 ~ 6) - /?] Q 



n = [Mfrf + hi (g + l)(6-5 7 ) + /?] {4^iV| [2( (7 - 1) + 7 2 ] + (7 - 1) P 2 } 



v 1 = - 



6/i 2 7 3 /3 



, v 2 = f 3 = C 1/2 7 1 



2/iiiV 3 7 

,4 + 5 

where we have set: 

A = 3/ii7|l + 3/ii||/ii(37- 2)- 57 + 6| x 

x{( 7 - 1) [h\ (7 - 1) (7 7 - 6) 2 + 2h! (7 - 2) (27 7 2 - 37 7 + 6) + 
+ ( 7 -2) 2 (9 7 -l)]} 1/2 

B = I8hj~/ 2 (7 - 1) (37 - 2) + 3/i? (66 7 4 - 427 7 3 + 8O87 2 - 588 7 + 144) - 
-3hj (1O67 4 - 595 7 3 + I2OO7 2 - 1052 7 + 336) - 

-3/ii (9O7 4 - 499 7 3 + 984 7 2 - 844 7 + 272) + (3 7 - 2) (35 7 - 36) (7 - 2)] 

A = 2[27/i 4 7 (7 - l) 2 (37 - 2) - I8/1? (157 4 - 62 7 3 + 937 2 - 587 + 12) + 
+3hj (757 4 - 384 7 3 + 704 7 2 - 56O7 + 168) - 
-6/ii (3O7 3 - 127 7 2 + I667 - 68) + (35 7 - 36) (7 - 2)] 

P = 2q - 3 7 + 2, 



10 



c = P 2 (7-2)(l + 3/ti)(g+l) 

4/niVf {3/if (4g 7 - /3) - 4/m [q (5 7 - 4) + 2 (7 - 2)] + 3/3} 

We note that the above solution is defined only when h± > — 1 and belongs to the subclass 
= 0. In addition the state parameter satisfies: 

1 u 1 1 2 (3-/ii) 

-1 < /ii < -- =>- 1 < 7 < 



3 ' 5 - 3/ii 

1 , 2(3-/ii) 3 

-3 </ii<0 ^T^r<^<2- 

It was not proved possible to study the stability properties of the solution in the full state space. 
However the simpler problem of the subclass of models satisfying = can be treated. In 
particular a preliminary analysis indicates that this equilibrium point has a five dimensional 
stable manifold and may act as the future attractor at least for the models satisfying = 0. 

8. Type Vlh Extreme tilted Set of Equilibrium Points C ex t re m,e(V Ih) (h = —h\ and h\ 7^ —1/3): 

N 2 = N 2 ,N 3 = N 3 , v a v a = l, 



2 _ 24}nV-N 2 N 3 - 12hjN 2 N 3 + {N 2 - N 3 ) 2 + 12 

12 



S 13 = S 3 = S_=0, £+ = -hW-N 2 N 3 -l, Si = ^- (N 3 - N 2 ) 



Vl 



1, 9 = 2(1 + ^^-^2^3), 



_ -12hW-N 2 N 3 + l2h\N 2 N 3 - (N 2 - N 3 ) 2 

" = , U < 7 < 2 



where h\ < 0. 



3.3 The exact form of the self-similar equilibrium points 

Apart from the cases with extreme tilt which are not representing by exact solutions of the 
EFE, the rest equilibrium points, whenever exist, correspond to self-similar models i.e. models 
admitting a proper Homothetic Vector Field (HVF) H acting simply transitively on space-time: 

£H9ab = tyQab (3.1) 
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where ip =const. is the homothetic factor which represents the (constant) scale transformation 
of the geometrical and dynamical variables. 

For vacuum and non-tilted models, the corresponding self-similar solutions are all known (see 
e.g. Tables 9.1. and 9.2., pages 187-188 of |3|) whereas for tilted models the known self-similar 
solutions are of type II JO] and VIo ^3 • Taking into account the results of the previous 
section we conclude that only in type VI^ (i.e. h ^ — 1) there exists a self-similar model. In 
order to find the exact form of this new family of solutions we follow the procedure that was 
used in |19j (a complete study of the self-similar SH models of class B is reported elsewhere 

Adopting the notation of I21| . the KVFs {A^ a } and the dual basis {uj a } are: 

X 1 = d y , X 2 = d z , X 3 = d x + yd y + h R zd z . 

u 1 = e~ x dy, J 2 = e~ hRX dz, to 3 = dx. 
where the group parameter h R is formally related with h via (hi < 0): 

1 — /ii , 
h R = ±— . (3.2) 

It follows from the non- vanishing structure constants of the isometry group = C 23 = 1 
and the Jacobi identities, that the remaining non- vanishing structure constants C^ 4 of the simply 
transitive homothety group are: 

! _ h R (2-pi) 2 _ , 

P2 ~ 2 

and the HVF is: 

H = -i<9 t H ydy + zd z . (3.3) 

P2 ~ 2 2 - p 2 

In addition using the fact that in every perfect fluid model, the fluid velocity u a is conformally 
mapped by the HVF i.e. £n u a = ^> u a in conjunction with the EFE, we write out explicitly the 
frame components of the four-velocity and the self-similar metric: 

Fluid velocity 

Ai = «it P1-1 , A 2 = 0, A 3 = v 3 t (3.4) 



Metric 

/ c n t 2 ^-^ c 13 t Pl \ 

g a p = C22t 2(h R - P2 +2)/h R (3 5) 

V ci 3 t pi c 33 t 2 J 

where u a = T (S^ + A a u>^) and g = —dt 2 + g aj 3U>^uj^dx a dx b . 

The various integration constants appearing in (|3.4j) and ()3.5|) are given by the following 
expressions: 

2 _2h 2 R (s- 1) + 2h R (2s - 1) +p 2 - 2 

7 "27TT' Pl ~ hi 
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ess = c n (h R + 2)(2s + l) [2h 2 R s + 2h R (2s - 1) + p 2 - 2] + 4//Q^iT 
ci 3 4cn/io/i|fir 



r 2 = cn(2a + 1) [h% (p 2 + 2s - 2) + 2/t fl (2a - 1) +p 2 - 2] 

2/Io4 (C11U3 - Cl 3 Vl) 



V3 



h R v\ [2h 2 R s - 2h R (pa - 2s - 1) + p 2 - 2] 
2c u s (h R + 2) [2/i|s + 2h R (2s - 1) + p 2 - 2] 



(2s + 1) {4/i|s (2s - 1) (2p 2 + 2s - 3) 



-24 



+ 



+/> 3 



+4 



+ 



(5s - 1) - p 2 (24s 2 - 3s - l) - 2s (20s 2 - 32s + 9) 
p 2 2 + 4p 2 (l2s 2 - 10s + l) + 4 (32s 3 - 48s 2 + 20s - l) 
p\ (16s - 5) + 12p 2 (4s 2 - 6s + l) + 4 (l6s 3 - 40s 2 + 24s - 1 
+2h R (pa - 2) [pa (2s + 1) + 2 (8s 2 - 6s - l)] + 4s (p 2 - 2) 2 } x 
x {24 (h R + 2) (s - 1) [2h R (pa - 1) - pa + 2]} _1 



+ 



C13 (fa + 2) [2/t fl (2s - 1) + p 2 - 2} [2h 2 R s + 2h R (2s - 1) + p 2 - 2] 
24 [4 (pa + 2s - 2) + 2h R (2s - 1) +p 2 - 2] 



c u = [2h R (2s - 1) + p 2 - 2] \ 2h R s + 2h R (2s - 1) + p 2 - 2 
x{44« (2s - 1) \p2 (s + 1) + 2s 2 - 2s - 1 
-2h\\p\ (3s 2 + l) - p 2 (lOs 3 + 18s 2 - 13s + 5) - 
-2 (24s 4 - 34s 3 + 4s 2 + 5s - 2)] - h R \2p\ (s - 1) + 
+p\ (8s 2 - 24s + 15) - 4p 2 (2s 3 + 18s 2 - 28s + ll) 



13 



_4 (48s 4 - 76s 3 + 24s 2 + 14s - 7)] - h 3 R [3p 3 2 (s - 1) + 
+Pl (4s 2 - 38s + 23) - Ap 2 (l2s 3 - 28s + 13) - 
-4 (32s 4 - 88s 3 + 60s 2 - 5)] + 2h\ (2 - p 2 ) x 
x [p 2 (6s 2 - 13s + 4) - 2 (24s 3 - 34s 2 + 13s - 2)] + 
+h R (p 2 - 2) 2 [2 (l2s 2 - 13s + 3) - 3p 2 (s - 1)] + 
2( P2 -2) 3 (s-l)} x 

x{8h R s (s - 1) [h 2 R (p 2 + 2s - 2) + 2h R (2s - 1) + p 2 - 2] 3 }" 1 

p 2 = 2{h R \(h R + 2)[h R (2s-l)-s]\x 

x [(2s - 1) (h 2 R (s + l) 2 (2s - 1) + 2h R (4s 3 - 4s 2 + 5s - l) + 

+8s 3 - 28s 2 + 10s - 1)] 1/2 + h\ (s + 1) (2s - l) 2 - 

-h% (2s 3 - 21s 2 + 17s - 4) - h R (8s 3 + 9s - 4) + 

+h 2 R (20s 3 - 16s 2 + 4s - l) - 2h R (4s 3 + 6s 2 - 8s + l) + 4s 2 } > 

x{2h R (8s 2 - 5s + l) - h 3 R (l6s 2 + 6s - 5) + 

+h 2 R (20s 2 - 14s + l) - 2h R (8s 2 - 8s + l) + 4s 2 }" 1 



From the expressions of the constants p 2 , T and the positivity of the energy density it can be 

verified that the above family of self-similar tilted perfect fluid solutions of type VI^ is defined 
for 

6 <s< 2i^T T^r <7< 2 (3 - 6) 

when — 1 < ha < — \ and 

hR 1 2 (3 -hi) , 

< a < ^ 1 < 7 < i ^ (3.7) 



2h R -l 2 5 - 3/ii 
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when — ^ < h R < 0. We note that the inequalities involving the bounds of the group parameter 
h R are reduced to — | < h\ < and — 1 < hi < — | respectively as expected from the analysis 
in the previous section. 

As a concrete and illustrative example of the general results given above, we present the 
self-similar metric and the associated energy density for the case of a radiation fluid i.e. 7 = 

3 w 5 4 . 



\(h R + 2) (2h R + 1)| y/25h? R -Ah R + A- I2h% - 58h% + 7h 2 R - 32h R - 4 
Pl = 4 (2h R + Wh R - 5h 2 R + 4h R + 1) 



+ 2) (2/ir + 1)| \J2hh R - 4h R + 4 + 10/i|, + 33/ifj + 52/i| - 22/i| + 36h R + 8 
P2 ~ 4h R 2 (2h% + lO/i^ - 5h? R + Ah R + l) 

/i = {-3fc B + 2| (6/t| - 5/i| + 6h 2 R + 6/ijj + 2) ^25h 2 R -Ah R + A + 

+30h% - 29/4 + 16/l l - 30/l l - 10/l fl + 2} x 

x{4t 2 (2/1^ + 10h% - 5h 2 R + 4h R + l) [-fo (2/i fi - 1) ^25h 2 R -4h R + 4 + 



+Wh% - 5h 2 R + 6h R + 2]} 



-1 



It turns out that the energy density is increasing as the group parameter varies in the interval 
h R € ^—1,-7^ or in terms of the "original" parameter h £ ^— §,0j. 

4 Concluding Remarks 

A long term goal of the qualitative study of SH tilted perfect fluid models rely on the determi- 
nation of the equilibrium points of their dynamical state space in order to be able to make solid 
conjectures regarding the dynamics, at the asymptotic regimes, of the corresponding models. 
In the present paper, by exploiting the orthonormal frame formalism and the general form of 
evolution equations given in t lT|, we have found the complete set of equations describing the 
evolutionary behaviour of the important class of Bianchi type VI/j tilted perfect fluid models. 

As a result we have identified the dynamical state space to be a bounded region T> C R 
subjected to the constraints (|2,24|) - (|2.26[) and we have given all the equilibrium points (both 
vacuum and non- vacuum) . This study led us to a new family of exact solutions of EFE which 
are represented by the self-similar metrics (j3.5|) . This family of rotational tilted perfect fluid 
models, belongs to the subclass = and has some interesting stability properties. 

In particular we have seen that this model is defined only for — 1 < h < and has a five 
dimensional stable manifold in the ranges h x < -± 1 < 7 < 2 g_~y and -§ < hi < =^ 

2 5-3fei' > ^ 7 ^ §■ ^ n con j unc tion with the fact that for — 1 < h < the Collins non-tilted solution 

^(Vh) is stable whenever § < 7 < 2 fc^ } OH we conclude that, at the value 7 — ^ 

2 We note that, in contrast with the case of VIo models in which a similar situation occurs concerning the value 
7 = I , there is no acceptable tilted solution for the value 7 = • 
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and when 1 + 2>h\ > 0, there is an exchange of stability between the Collins model and the 
present solution. Moreover for 1 + 3h\ < both models are stable whenever 1 < 7 < "5^3^ 

and they are also destabilised at 7 = 2 ^~^ . However, due to the complexity of the situation, 
there are no straightforward conclusions regarding the asymptotic behaviour of the models in 
the full state space. Perhaps a more sophisticated choice of the state variables is needed in order 
to efficiently answer many open questions. Nevertheless the case h = — 1 appears to be more 
tractable, permitting a complete analysis of its dynamical properties. We also pointed out that 
the above procedure can be used to determine the equilibrium points for the rest of the Bianchi 
models and for the case of the exceptional models h = — g as well. These matters will be the 
subject of a future work. 

A quick glance of the equilibrium points in section 3 shows the lack of physically acceptable 
models for h = — 1 which proves (this fact can be confirmed using the geometric results of |20| ^ 
the non-existence of tilted self-similar type III models. Furthermore in all (non extreme tilted) 
cases, the tilted fluid velocity is not parallel to A a . According to Theorem 3.2. of [2] this implies 
that self-similar type VI^ models have necessary non-zero vorticity. In conclusion we have the 
following proposition: 

Proposition 1. There are no type III and irrotational type Vlh tilted perfect fluid models 
which admit a proper HVF. 

We finally note that a stability analysis of the irrotational tilted perfect fluid models of type 
Vlh shows that the equilibrium points (both non-tilted and tilted) found in section 3 are future 
unstable which, by means of Proposition 1, implies that irrotational type Vlh models are not 
asymptotically self-similar but rather they are asymptotic to an extreme tilted model. 
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